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We compute the 0{a^s) QCD corrections to the fully differential cross-section pp WX IvX, 
retaining all effects from spin correlations. The knowledge of these corrections makes it possible to 
calculate with high precision the W boson production rate and acceptance at the LHC, subject to 
realistic cuts on the lepton and missing energy distributions. For certain choices of cuts we find large 
corrections when going from next-to-leading order (NLO) to next-to-next-to-leading order (NNLO) 
in perturbation theory. These corrections are significantly larger than those obtained by parton- 
shower event generators merged with NLO calculations. Our calculation may be used to assess and 
significantly reduce the QCD uncertainties in the many studies of W boson production planned at 
the LHC. 
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Production of clectroweak gauge bosons is a vital com- 
ponent of the hadron collider physics program. The large 
production rates for this channel at the LHC will facil- 
itate several important precision measurements. LHC 
experiments plan to determine the W boson mass and 
width with errors of 5Mw ~ 15 MeV and 5Tw ~ 
50 MeV, respectively Q. The Weinberg angle svclOw 
can be extracted from the forward-backward asymmetry 
of the lepton pair m. pp ^ Z ^ e+e^ with a precision of 
1 X 10~^. The precision possible in these channels at high 
luminosities makes these measurements competitive with 
LEP results. Searching for deviations from predictions 
in di-lepton events with large invariant mass, missing en- 
ergy, or transverse momentum probes extensions of the 
Standard Model which contain new gauge bosons. 

In addition to its high rate, electroweak gauge boson 
production has a simple, distinct experimental signature. 
This also makes it a useful process for calibrating and 
monitoring machine and detector performance. Z and 
W production can be used to determine and monitor the 
hadronic and partonic luminosities at the LHC . This 
requires a theoretical prediction for the cross section to 
the highest possible precision, since this error propagates 
into all other measurements through the luminosity un- 
certainty. Determination of the LHC luminosity to 1% 
accuracy is the ultimate goal of this procedure 0] ■ This 
sets the precision required for theoretical predictions. 

When the desired precision on the production cross 



section is at the few percent level, many subtle effects 
must be included. Both 0{a) electroweak effects and 
©(a^) QCD effects must be calculated. The electroweak 
corrections to pp ^ W ^ Iv were computed in |^ , where 
the importance of final state photon radiation in the W 
decay was observed. NLO computations of the QCD cor- 
rections to electroweak gauge boson production were first 
obtained in the late seventies Q. The W boson momen- 
tum distribution was investigated via resummation tech- 
niques in [s^. Currently, the NNLO QCD corrections are 
known for both the inclusive production cross section |^ 
and for the gauge boson rapidity distributions T\. The 
NNLO corrections are typically in the few percent range 
at the LHC, and must be included in both precision elec- 
troweak studies and the luminosity determination. 

Existing calculations of the NNLO QCD corrections 
to this process do not include all effects needed for a 
percent level theoretical prediction. Phenomenological 
applications of Z and W production require significant 
cuts on the phase-space of the final state leptons. For ex- 
ample, all LHC experiments will impose constraints on 
the transverse momenta and rapidities of the final state 
charged leptons. Cuts on the missing energy will also be 
employed to identify the neutrino from W decay. Cal- 
culations that treat the leptons inclusively are therefore 
not fully realistic. They can be used to make estimates, 
but they are not sufficient for precision measurements. 

The calculation of the full NNLO QCD corrections is 
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complicated by the spin-one nature of the gauge bosons. 
If they were spin-zero bosons, fully differential results 
could be obtained from Ref. 0, where the rapidity dis- 
tributions for Z and W bosons were computed through 
NNLO. That result could be combined with the known 
double differential distribution in transverse momentum 
and rapidity to fully determine the gauge boson kine- 
matics. The decay of spin-zero bosons in their rest frame 
is isotropic, and the final state distribution of interest 
could be obtained by assuming a flat decay distribution 
in the gauge boson rest frame. There would be no cor- 
relation between the production and decay of the boson. 
However, since the Z and W are spin-one bosons, the 
Z l^l^ and W ^ Iv decays are not isotropic, and 
there are " spin correlations" between the production and 
decay channels. The importance of spin correlations for 
the production of gauge bosons at the LHC was recently 
emphasized in f9| . 

The computation of the differential cross section for 
pp -> WX -> IvX through NNLO in QCD is a difh- 
cult theoretical challenge. While techniques for perform- 
ing differential NLO calculations have been known for 
many years the corresponding technology for ob- 
taining NNLO results is still in its infancy. In a recent 
series of papers |l2j| , we have developed a method for per- 
forming these calculations. This technique features an 
automated extraction of infrared singularities from the 
real radiation matrix elements and a numerical cancel- 
lation of these divergences with the virtual corrections. 
We describe below the application of this method to the 
computation of the W production cross section with all 
spin correlations included. 

We compute the partonic cross sections ij — > IvX as 
perturbative expansions in the strong coupling constant 
as- We specialize here to W~ production. At LO, the 
W~ boson is produced in the collision of an up-type anti- 
quark and a down-type quark. At NLO, gluon-quark and 
gluon-antiquark scatterings also contribute. A variety of 
partonic processes contribute at NNLO. These have been 
enumerated in great detail in . In our discussion below 
we use ud — > W~ X — > evX as an example, since it con- 
tains all the complexities present in the full calculation. 
All partonic channels have been included in our result. 

There are three distinct contributions contained in the 
ud initiated process: the two-loop virtual corrections, the 
one-loop virtual corrections to single gluon emission, and 
tree-level double-real radiation processes with two addi- 
tional partons in the final state. These must be combined 
in the presence of an infrared-safe measurement function 
to produce a finite result. We use dimensional regular- 
ization to regulate all ultraviolet, soft, and coUinear di- 
vergences. 

The two-loop virtual corrections to the ud process are 
straightforward to compute. They are very similar to 
the 0{a^g) corrections to the quark form factor studied 
in pTl |: however, the W production calculation must in- 
clude the two-loop corrections to non-singlet axial cur- 
rent. Care must be taken to define this correctly in 



d = 4— 2e dimensions; we discuss this issue further below. 
We use the implementation of the Laporta algorithm 
described in Ref. 0| to reduce all required two-loop in- 
tegrals to a minimal set of master integrals. The master 
integrals needed for this computation are well known . 

We obtain the one-loop correction to the single gluon 
emission process ud — > eueg using a combination of two 
methods. We first use the Laporta algorithm to express 
all one-loop Fcynman integrals relevant for this process 
through master integrals. The master integrals must be 
integrated over the final state phase-space subject to the 
kinematic constraints under consideration. It is not pos- 
sible to perform this integration analytically, since we 
want an expression valid for arbitrary cuts. Numerical 
integration is also not straightforward because of soft 
and coUinear singularities. We employ the method devel- 
oped in |l2j| to extract the singularities in a constraint- 
independent way as poles in e before integrating over the 
phase-space numerically. This technique maps the final 
state phase-space onto the unit hypercube and uses it- 
erated sector decomposition ITsI to extract all soft and 
coUinear singularities. 

We use essentially the same algorithm to compute the 
double-real radiation corrections. A detailed description 
of this method, which studies in detail both the one and 
two parton emission corrections, can be found in [T^ . 

We now discuss a few new features of this calcula- 
tion, first explaining how we treat the axial current in d- 
dimcnsions. This issue arises from Dirac structures of the 
form TrH[r(i)75]TrL[r(2)75], where T^^'^) denote generic 
products of Dirac matrices and TrH,L refer to traces over 
hadronic and leptonic degrees of freedom, respectively. 
These traces do not vanish when the final state phase- 
space is sufficiently constrained. A consistent extension 
of the axial current to d-dimensions is given in T^. It 
utilizes an anti-commuting 75 and contains additional 
renormalizations relative to the vector current in order 
to maintain the Ward identities. We use this prescrip- 
tion in our calculation. 

Even after all three components of the hard-scattering 
cross section are combined, coUinear counterterms are 
needed to remove initial state coUinear singularities. 
In these coUinear counterterms were treated analyt- 
ically. Such an approach is not sufficiently flexible to 
handle cuts on the W decay products. However, it is 
straightforward to extend the numerical approach used 
for the other NNLO components in to obtain the 
desired results. 

We have essentially two checks on our calculation. 
First, considering different cuts on the electron trans- 
verse momentum and rapidity as well as on the missing 
energy, we verify cancellation of the divergences in the 
W~ production cross section. Because the divergences 
start at 1/e^ at NNLO, the cancellation of all divergences 
through 1/e provides a stringent check on the calculation. 
We also check that the vector and axial contributions are 
separately finite, as required. A second check is obtained 
by integrating fully over the final state phase-space and 



comparing against known results for the inclusive cross 
section. We find excellent agreement with the results 
of for all partonic channels. 

We now discuss the results of our calculation. We first 
present the input parameters. We use the MRST parton 
distribution functions |17| at the appropriate order in as- 
We use mw = 80.451 GeV and work in the narrow width 
approximation, although this restriction can be easily re- 
moved. We set iKdl = 0.974, = \Vcd\ = 0.219, and 
|Vcs| = 0.996, and obtain \Vab\ and \Vch\ from unitar- 
ity of the CKM matrix. We neglect contributions from 
the top quark; these have been shown to be small in the 
inclusive cross section 1^ . For electroweak input param- 
eters, we use sin^ 9w — 0.2216, Q!qed('7iz) = 1/128, and 
Br(Vl^ — > ev) = 0.1068. We set the factorization and 
renormalization scales to a common value, = = 
and employ various choices of /i in our numerical study. 

We find that NNLO corrections depend on the cuts and 
can change rapidly from very small to fairly substantial. 
We consider cuts of the form 

p1>p];""", 177-= I < 2.5, Sf'^'^ > 20 GeV, (1) 

and use the values p"^""" = 20,30,40,50 GeV. The 
choices p^'"'" = 20 and 40 GeV were considered in the 
study of (9j. = 20 GeV is similar to cuts that 

will be employed by the ATLAS and CMS collaborations, 
while p^""" = 40 GeV was chosen in to illustrate the 
potential sensitivity of the QCD radiative corrections to 
experimental cuts. We first present results for the lepton 
invariant mass distribution for on-shell W~ production 
in TableQI to give a feeling for the magnitude of the cross 
section for each set of cuts. The numerical precision for 
all NNLO numbers is 1% or better. We note that the row 
labeled "Inc" denotes the fully inclusive cross section. 



e.min 

Pi 


LO 


NLO 


NNLO 


Inc 


11.70,13.74,15.65 


16.31,16.82,17.30 


16.31, 16.40, 16.50 


20 


5.85,6.96,8.01 


7.94,8.21,8.46 


8.10,8.07,8.10 


30 


4.305, 5.12,5.89 


6.18,6.36,6.54 


6.18,6.17,6.22 


40 


0.628,0.746,0.859 


2.07,2.10,2.11 


2.62,2.54,2.50 


50 


0,0,0 


0.509,0.497,0.480 


0.697,0.651,0.639 



TABLE L The lepton invariant mass distribution da/dM^ , 
M — mw, for on-shell W production in the reaction pp 
W-X e-vW, in pb/GeV^ for various choices of p"^'™'", 
GeV and = mw /'2, mw, 2mw- 

There are a few things to notice about these num- 
bers. First, at LO, there is no additional hadronic ra- 
diation in the final state for the lepton and neutrino to 
recoil against, so the transverse momentum is restricted 
to < TOH//2. The cross section is therefore very small 
for p^""" = 40 GeV, and vanishes for p*]^""" = 50 GeV. 
This restriction is lifted at NLO when there is an addi- 
tional parton for the W to recoil against. Very near this 
boundary, the width of the W can be an important ef- 
fect. It will induce a (tiny) cross section for p*]^™'" = 50 
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GeV at LO, and it will shift the result for p^™'" = 40 
GeV since this value is close to mn//2. We have not in- 
cluded the M^-width in our results. However, we have 
checked using the results in that finite width effects 
change the acceptance by only 7% at LO, and by less at 
NLO. We are therefore confident that our discussion and 
conclusions are not affected by this omission. 

Another feature to notice is that the corrections are 
large for higher choices of p']^'"'" and that the dependence 
on p^™'" is strong. For example, for the scale choice 
/i = mw, we observe a 22% increase when going from LO 
to NLO inclusive cross section, followed by a decrease of 
2.5% when NNLO corrections are included. We obtain 
similar results for p^""" = 20, 30 GcV. However, the pat- 
tern of corrections is much different for p*^™'" — 40, 50 
GeV; we find corrections of 18-27% for p^j;""" = 40 GeV 
and 33-37% for 50 GcV when going from NLO to NNLO, 
depending on the choice of scale. 

The remaining scale dependences can be seen from 
Table We define the scale dependences of 

the cross section ax as Actx = 2(Max[(Tx(Ai)] — 
Min[crx(M)])/(Max[crx(M)]+Min[crx(M)]), where we take 
the maximum and minimum values from among the three 
studied scale choices. Aax therefore gives the scale vari- 
ation uncertainty band for the observable ax ■ The scale 
dependence is reduced to the percent level or less for 
the inclusive case, and for the cuts p^™'" = 20, 30 GeV. 
Moreover, the NNLO results lie within the NLO uncer- 
tainty bands. This is not the case for the other two 
choices of p^"^™; here, the scale dependence actually in- 
creases to the 5-8% level at NNLO, and the NLO scale 
dependence completely underestimates the higher-order 
radiative corrections. This is not completely unexpected, 
since for these values additional partons to recoil against 
only appear at NLO. The NNLO results therefore serve 
as the first radiative corrections for these p^™'" choices. 
However, it indicates the care that must be taken when 
using the scale variation as a measure of the theoretical 
error. 





y4(NL0) 


A(NNLO) 


20 


0.487,0.488,0.489 


0.497,0.492,0.491 


30 


0.379,0.378,0.378 


0.379,0.376,0.377 


40 


0.127,0.125,0.122 


0.161,0.155,0.152 


50 


0.0312,0.0295,0.0277 


0.0427,0.0397,0.0387 



TABLE II: Acceptances at NLO and NNLO for various 
choices of p*!^™'" and jj, = mw /2,mw j'imw ■ 



Another important quantity to study is the experimen- 
tal acceptance, defined as the ratio of the cross section 
after cuts over the inclusive cross section. We present the 
acceptances at NLO and NNLO in Tabic ^1 We again 
note that for the choices p^™'" — 40, 50 GeV, the NLO 
scale dependences completely underestimate the NNLO 
corrections. The NNLO shifts in acceptances are very 
large for p'j;™'' = 40, 50 GeV, reaching 25% for 40 GeV 
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and 40% for 50 GeV. For the other choices of cuts, the 
NNLO acceptances are identical to the NLO ones within 
numerical errors, indicating stabilization of the pertur- 
bative expansion. 

The cross sections and acceptances for the transverse 
momentum cuts p'j^""" = 20, 40 GeV were recently stud- 
ied in 9]. The primary tool used in that analysis was 
the Monte Carlo event generator MC@NLO, which con- 
sistently combines NLO corrections with the HERWIG 
parton shower . In , MC@NLO is used to estimate 
the importance of QCD effects beyond those included in 
NLO calculations. Differences between the cross sections 
and acceptances of a few percent are found when com- 
paring NLO and MC@NLO for both p^^'"'" choices. The 
authors of then conclude that higher order corrections 
beyond those in MC@NLO are unlikely to change signif- 
icantly the results they find. 

We believe that the few percent differences between 
MC@NLO and NLO cannot be used as an estimate of 
higher order corrections. This is because few percent 
shifts coming from hard emissions generically occur in 
process-dependent radiative corrections that cannot be 
described by parton showers. We can support this asser- 
tion with the following observations. It follows from Ta- 
ble 2 in |§| that adding the parton shower to the LO cross 
section fails to properly predict the NLO cross section. 
For p^""" — 20 GeV, adding the parton shower to the LO 
result decreases the cross section by 8%, while the NLO 
correction increases it by 4%. While both HERWIG and 
NLO corrections increase the LO result for — 40 

GeV, the magnitude of the shifts differ by 45% relative to 
the LO cross section. In addition, the NNLO corrections 
to the cross-sections presented in Table differ signifi- 



cantly from the estimate of these corrections in p| . We 
conclude that the results for the cross sections and ac- 
ceptances obtained in cannot be used if few percent 
precision is required. This is particularly true for cuts 
where hard gluon emissions are expected to be large, such 
as for p'j]™'" = 40, 50 GeV. We note that since the correc- 
tions are large for p*]^""" — 50 GeV, which is well above 
the LO kinematic boundary for the electron transverse 
momentum at mvi//2, our results are not caused by large 
logarithms which spoil the perturbative expansion. 

In this Letter we report on the computation of the 
NNLO QCD corrections to the fully differential cross sec- 
tion pp WX IvX at the LHC, with all spin cor- 
relations included. For inclusive enough cuts, this cal- 
culation provides the percent-level theoretical accuracy 
needed for the W cross section when realistic experimen- 
tal cuts are imposed on the final state leptons. We find 
that the QCD corrections exhibit significant dependence 
on the lepton minimum transverse momentum. For high 
values of this cut, the corrections may be very different 
than the inclusive NNLO results. Our calculation can 
be easily extended to include Z production, finite width 
effects, and pp collisions. The last case is particularly in- 
teresting because of its importance for the Tevaton Run 
II physics program. These extensions will be discussed 
in detail in a forthcoming publication. 
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